SPECIFICATION 

METHOD FOR MEASURING DIGGING POSITION 



BACKGROUND OF THE INVENTION 
[Technical Field of the Invention] 

The present invention relates to a method for determining the digging 
position in a non-open-cut method of excavation and, more particularly, to a 
method which ensures accuracy in detemiining positions by lessening die 
influence of a noise magnetic field of frequency components close to that of die 
signal magnetic field to be measured. 

[Prior Art] 

A horizontal drilling mefliod, which is one of the non-open-cut methods of 
this kind, uses a small-diameter pipe of 100 mm or less across for horizontally 
digging in die ground, and accordingly, such a precision position determining 
apparatus as used in an ordinary small-diameter driving metiiod of excavation 
cannot be placed near a drill. To solve tiiis problem, it is customary in die art to 
adopt a method in which an AC magnetic field is generated by a coil mounted in 
die drill head and detected by an above-ground magnetic sensor like a coil to 
determine die current digging position. 

This method is simple and easy, but since die magnetic field by the coil is 
a dipole magnetic field, it rapidly attenuates widi distance from die cod. Hence, 
diis method has a defect of inability to achieve liigh-reliability determination of die 
digging position when a power line or simdar magnetic noise source is present in 
die vicinity of die place where to perfomi the position detennination. 

SUMMAR OF THE INVENTION 
An object of die present invention is to provide a position determination 
method diat permits high-reliability detennination of die digging position by 
detecting on the ground the AC signal magnetic, field provided from a cod housed 
in die drill head even if a noise magnetic field is present which affects die position 
detennination. 

To attain die above object, a digging position detemiining metiiod 
according to die present invention for non-open-cut excavation, which senses an 
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AC magnetic field provided from a magnetic field source by a magnetic sensor 
provided on the ground and calculates the position of die magnetic field source 
from the magnitude and direction of die sensed magnetic field, said method having 
a construction characterized in that 
5 In case where in addition to a signal magnetic field generated by said 

magnetic field source, there exists a noise magnetic field generated by a nearby 
current, 

at lease one of the position of said magnetic field source, die tilt angle of 
said magnetic field source to the vertical direction and the azimuth angle of its 

10 axial direction of said magnetic field source in a horizontal plane is calculated, 
from a projective component of die magnetic field sensed by said magnetic sensor 
and projected on a plane or straight line orthogonal to a vector-valued direction of 
said noise magnetic field. 

That is, in 1999 year's investigation and research relating to useful 

15 utilization techniques of energy resources: entitled cc Research for low-loss 
energizing techniques in establishment of advanced telecommunication 
network",(executed by Composite Development System for New Energy 
Industrial Technique), it is described that external noise magnetic fields, which 
affects the position determination in the non-open-cut method of excavation, is 

20 mostly generated by a current of some kind hi this case, the magnitude of die 
noise magnetic field varies irregularly witii time, but its vector-valued direction is 
constant at each field sensing position. 

The present invention attains its object dirough die adoption of the 
following steps (A) and (B). 
25 (A) The direction of the noise magnetic field is detected, and a sensed 

magnetic field in which die noise magnetic field and a signal magnetic field are 
mixed is projected on a plane or straight line orthogonal to the direction of die 
noise magnetic field to obtain a projective component. 

Since die projective component is tiieoretically free from a component 
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derived from the noise magnetic field, at lease one of the position, azimuth angle 
and tilt angle of the magnetic field source is calculated so that the magnitude of the 
projective component (in the case of projection on the straight line) or it magnitude 
and direction (in the case of projection on the plane) is substantially equal to a 
theoretically calculated value of a corresponding quantity of a magnetic field 
generated from a magnetic field source, or ,has a minimum difference between the 
former and the latter. The sensed magnetic field is obtained by sensing magnetic 
fields at different positions whose number is determined by how many ones of the 
position, azimuth angle and tilt angle of the magnetic field source are unknowns 
and how such unknowns are calculated. 

(B) To obtain die noise magnetic field, 

a) When the number of noise magnetic fields is virtually one: 

• The noise magnetic field is sensed to obtain its direction essentially in the 
absence of a signal magnetic field. 

• When die noise magnetic field has frequency components at frequencies 
different from that of the signal magnetic field, the frequency components are 
measured to obtain the direction of the noise magnetic field. 

b) When the number of magnetic fields is virtually two: 

The frequency components of a first one of the two noise magnetic fields, 
which are widely spaced from the frequency components of the second noise 
magnetic field and Hie signal magnetic field, are measured to obtain a vector 
direction of the first noise magnetic field, and the frequency components of die 
second noise magnetic field, which are widely spaced from the frequency 
components of the fust noise magnetic field and the signal magnetic field, are 
measured to obtain a vector direction of the second noise magnetic field. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The features of die present invention will be clearly understood from the 
following description taken in conjunction with the accompanying drawings, in 
which: 



Fig. I is a perspective view explanatory of the placement of a magnetic 
sensor in the present invention method; 

Fig. 2 is a vector diagram explanatory of the principle of measurement 
according to the present invention method in a case of one noise magnetic field 
5 source; 

Fig. 3 is a vector diagram explanatory of the principle of measurement 
according to the present invention method in a case of two noise magnetic field 
sources; 

Fig. 4 is a flowchart illustrating a procedure for obtaining a projective 
10 component according to the present invention method in a case of one noise 
magnetic field source; 

Fig. 5 is a flowchart illustrating a procedure for obtaining a projective 
component according to the present invention method in a case of two noise 
magnetic field sources; 

1 5 Fig. 6 is a flowchart illustrating a procedure for obtaining a signal magnetic 

field in a case of one noise magnetic field; 

Fig. 7 is a flowchart illustrating a procedure for obtaining a signal magnetic 
field in a case of two noise magnetic fields; 

Fig. 8 is a flowchart showing a measurement procedure in the present 
20 invention method; 

Fig. 9 is a flowchart showing a procedure for calculating the position of the 
signal magnetic field according to die present invention method in a case of one 
noise magnetic field when die number of unknowns and the number of equations 
are die same as to each oilier; 
25 Fig- 10 is a flowchart showing a procedure for calculating the position of 

die signal magnetic field source according to die present invention metiiod in a 
case of one noise magnetic field source when die number of equations is larger 
than the number of unknowns; 

Fig. 11 is a flowchart showing a procedure for calculating the position of 
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the signal magnetic field according to the present invention method in a case of 
two noise magnetic field sources when the number of unknowns and the number 
of equations are the same as to each other; 

Fig. 12 is a flowchart showing the procedure for calculating the position of 
5 die signal magnetic field source according to die present invention method in a 
case of two noise magnetic field sources when the number of equations is larger 
than the number of unknowns; 

Fig. 13 is a flowchart showing an example of a procedure for obtaining the 
signal magnetic field through the use of the magnitude of the signal magnetic field 
1 0 vector; 

Fig. 14 is a flowchart showing another example of a procedure for 
obtaining die signal magnetic field dirougli the use of die magnitude of the signal 
magnetic field vector; 

Fig. 15 is a perspective view explanatory of the placement of a magnetic 
15 sensor for determining the direction of a noise magnetic field according to the 

present invention metiiod when a digging head, which is a signal magnetic field 

source, is distant from the position of measurement; 

Fig. 16 is a flowchart showing a procedure for calculating the direction of 

die noise magnetic field according to die present invention method in a case where 
20 die digging head is distance from die position of measurement and die signal 

magnetic field source stops to generate die magnetic field; 

Fig. 17 is a flowchart showing a procedure for calculating the direction of 

die noise magnetic field according to die present invention mediod when the signal 

magnetic field and die noise magnetic field are mixed to each otiier; 
25 Fi & 18 is a signal frequency spectrum diagram explanatory of how the 

operation of selecting a frequency for maximizing a frequency spectrum is used in 

the process for calculating the direction of die noise magnetic field in die flowchart 

shown in Fig. 13; 

Fig. 19 is a flowchart explanatory of a first method by which the flow of 
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candidate vector calculating process is used in the process for calculating the 
direction of the noise magnetic field in the flowchart shown in Fig. 17; 

Fig. 20 is a flowchart explanatory of a second method by which die flow of 
candidate vector calculating process is used in the process for calculating die 
direction of the noise magnetic field in the flowchart shown in Fig. 17; 

Fig. 21 is a flowchart explanatory of a method by which the flow of 
process for evaluating the candidate vector and obtaining the direction of die noise 
magnetic field is used in the process for calculating the direction of the noise 
magnetic field in the flowchart shown in Fig. 17; 

Fig. 22 is a signal frequency spectrum diagram explanatory of how die 
operation of selecting a frequency for maximizing a frequency spectrum is used in 
die process for calculating the direction of the noise magnetic field in the flowchart 
shown in Fig. 13. 

Fig. 23 is a signal waveform diagram explanatory of the operation of 
specifying a period in which only a noise magnetic field exists tiirough utilization 
of a fact diat the amplitude of die sensed magnetic field signal becomes small 
during die OFF period of the signal magnetic field that is turned OFF by a 
predetennined procedure, in the calculation of die direction of the noise magnetic 
field according to die present invention mediod; 

Fig. 24 is a signal waveform diagram showing instantaneous variations in 
the amplitude of die sensed magnetic field signal when the signal magnetic field is 
periodically turned OFF in the calculation of the direction of the noise magnetic 
field according to die present invention method; 

Fig. 25 is a signal waveform diagram showing instantaneous variations in 
the amplitude of die sensed magnetic field signal when the signal magnetic field is 
randomly turned OFF in the calculation of the direction of die noise magnetic field 
according to the present invention method; 

Fig. 26 is a flowchart showing a method for calculating die direction of the 
noise magnetic field according to the present invention in which the signal 



magnetic field is turned OFF on the basis of a predetermined sequence, a period 
during which a particular magnetic field is OFF is specified by a sequence starting 
at the time when a correlation function between the sequence and the sensed 
magnetic field becomes maximum, and the direction of the sensed magnetic field 
5 in the specified period is regarded as the direction of the noise magnetic field; 

Fig. 27 is a flowchart showing a method for calculating the direction of the 
noise magnetic field according to die present invention in which the signal 
magnetic field is turned OFF on the basis of a predetermined sequence, and the 
stalling time of a sequence indicating the most likely ON/OFF state of the signal 
10 magnetic field is calculated from a plurality of times when a correlation function 
between the predetermined sequence and the sensed magnetic field becomes 
maximum; 

Fig. 28 is a flowchart showing the method for calculating he direction of 
the noise magnetic field according to the present invention in which the signal 

15 magnetic field is turned OFF on the basis of a predetermined sequence, the 
correlation function between the predetermined sequence and the sensed magnetic 
field is calculated at each of a plurality of time points at which Hie period of the 
sequence is equally divided, the magnetic field is projected at each time point to a 
vector formed by the calculated function, and one of vectors whose variance is 

20 minimum is regarded as the direction of the noise magnetic field; 

Fig. 29 is a flowchart showing another method for calculating the direction 
of the noise magnetic field according to the present invention in which the signal 
magnetic field is turned OFF on the basis of a predetermined sequence, a 
correlation function between the predetermined sequence and die sensed magnetic 

25 field is calculated at each of a plurality of time points at which die period of the 
sequence is equally divided, the magnetic field is projected at each time point to a 
vector formed by the calculated function, and one of vectors whose variance is 
minimum is regarded as die direction of the noise magnetic field; 

Fig. 30 is a flowchart showing die method for calculating die direction of 



the noise magnetic field according to the present invention in which the signal 
magnetic field is turned OFF on the basis of a predetermined sequence, and die 
starting time of a sequence indicating the most likely ON/OFF state of die signal 
magnetic field is calculated from a plurality of time points when the correlation 
5 function between die predetermined sequence and the sensed magnetic field 
becomes maximum; 

Fig. 31 is a flowchart showing still another method for calculating the 
direction of die noise magnetic field according to the present invention in which 
the signal magnetic field is turned OFF on the basis of a predetermined sequence, a 
1 0 correlation function between die predetermined sequence and the sensed magnetic 
field is calculated at each of a plurality of time points into which die period of the 
sequence is equally divided, die magnetic field is projected at each time point to a 
vector fomied by the calculated function, and one of vectors whose variance is 
minimum is regarded as the direction of the noise magnetic field; and 
15 Fig- 32 is a perspective view depicting an example of a magnetic field 

sensing frame for use in die present invention. 

DETAILED DESCRIPTION OF THE INVENTION 
As depicted in Fig. 1, in a case where the position of a digging head 2 
under die ground surface 1, which is a signal magnetic field source to be sensed, is 
20 determined using a magnetic sensor 4 placed on the ground surface 1 at a proper 
position when a power line or similar magnetic noise source 3, which generates a 
noise magnetic field, is placed near die digging position to be determined, there are 
present a signal magnetic field vector H s provided from die digging head 2 and a 
noise magnetic field vector H n from die magnetic noise source 3, such as power 
25 line. In diis case, die magnetic sensor 4 senses a vector H m tiiat is a combined 
version of die signal magnetic field vector H. and die noise magnetic field vector 
H„. 

Now, let die noise magnetic field of a position vector r and at time t be 
identified as a vector H n (r, t). On the odier hand, let die signal magnetic field 
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generated by magnetic filed generating means for position sensing be identified as 
a vector H s (r-r c , 0 C , t). Here, the vector 0 C is an angle of orientation of the 
magnetic field generating means, which is defined by three angles of rotation in the 
coordinate system fixed to the ground that is the coordinate system fixed to the 
magnetic field generating means. 

Since the noise magnetic field and the signal magnetic field are both 
sensed simultaneously and since the noise magnetic field vector H n (t, t) varies 
randomly with time, it is impossible to extract only the signal magnetic field vector 
H s (r-r c , 0 C , t) from the measured magnetic field vector H m (r-r e , 9 C , t) unless 
statistical properties of the noise magnetic field are known and the noise magnetic 
field is signal-wise orthogonal to the signal magnetic field. Even if die statistical 
properties for separating the noise magnetic field from the signal magnetic field are 
known prior to the determination of the position of die latter, the separation calls 
for a large amount of data, and hence die conventional scheme is not ever practical. 
According to die present invention, die direction of a vector e n (r) of die noise 
magnetic field vector H n (r, t) is obtained by separated means, and in a coordinate 
system shown in Fig. 2, a component vector H m p (r-r c , 9 C , t), projected to a plane 
vertical to die direction of a vector e^r) of die measured magnetic field vector 
H m (r-r c , 9 C , t) as shown in Fig. 4 (S 1 , S2, S3). 

HJXr-iv, 9 E , 0=H ra (r-r e) 0-(H m (r-r c> ')e n (r)> n (r). (1) 
This component does not contain the noise magnetic field for die reason given 
below. Since 

H,„ P (r-.;, *„. 0=H.(r- r c> 0 e , ,) + H n (r, 0=H,(r-r c , f ) + |H n (r, ^(r). (2) 
it follows diat 

25 H m P (r-r c , 0 C , ,)= H a (r- r c , 0 C> t)-^-^ 0-e n (r)>„(r). ( 3 ) 
from which it can be seen diat die projective component vector H m p (r-r c , 9 C , t) does 
not contain die component of the noise magnetic field vector H n (r, t). 

However, die projective component vector H ra p (r-r c , G c , t) has lost 
information of one axis by die projection on a plane vertical to die direction of the 
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vector e n (r). That is, since the same projective components are obtained 
irrespective of the magnitude of a component parallel to the vector e n (r), two 
independent components are obtained. 

Although any method can be used to obtain the two independent 
5 components; it is possible to use such a method as described below. 

That one of coordinate axes of a measurement coordinate system Cm 
(which will be described later on), which is not parallel to the direction e n (r) of the 
noise magnetic field vector H„(r, t), is chosen. Let a unit vector in die direction of 
the chosen coordinate axis by identified as a vector e m . A vector product, 

10 e p> ,=e m xe n (r), of die unit vector and the direction e„(r) is perpendicular to die 
direction of the vector e n (r), and hence it is contained in die plane of projection and 
is perpendicular to die coordinate axis e m . Let the magnitude of a vector obtained 
by projecting die projective component vector H m p (r-r c , 8 C , t) in die direction of die 
vector epj, including the direction of die vector, be represented (S4) by a value of 

15 HngV-fc, 9c, t). That is, 



The next step is to calculate a vector perpendicular to die directions of 
die vectors e^ and e n (r). The direction of the vector ep^ is also perpendicular to 

20 die direction of die vector e„(r), and hence it is contained in die plane of projection 
and is perpendicular to die direction of die vector e^i as well. Letting die 
projection of the projective component vector H m p (r-r c , 9 C , t) in tiiis direction be 
represented by H^Vr* 0 C , t), values H^Vr* 0 C , t)^ and H^Vr* 9 C , t)ep >2 
are two independent vectors into which the projective component vector H m p (r-r c , 

25 9 C , t) is separated. Here, 

(r - r c , ft, /)= HJV r c , 6> , /)- e p , 2 . (6) 
e P .2=e p . I xe n (r) (7) 

Then, by setting die position vector r c and angle-of-orientation 9 C of the 
magnetic field source so diat a theoretically calculated magnetic field, H^r-rc, 9 C , t) 
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generated by the magnetic field source of Ihe position vector r_substantially 
matches with a projective component H e p (r-r c , 8 C , t) on the same plane as that of 
die measured magnetic field, it is possible to detect the position and orientation of 
the magnetic field source. 

The above description has been given of a case where die number of noise 
magnetic fields is virtually one, but when die number of noise magnetic fields is 
two, letting die direction of two noise magnetic fields H nj (r, t), where j = 1, 2, be 
represented by vectors e,„(r) and e ll2 (r), use is made of die projection of die 
measured magnetic field on a direction vector eN(r)=e nl (r)x en2 (r) in the coordinate 
system of Fig. 3 as depicted in Fig. 5. That is, 

H,„ P (r-r c , 9 0 , ()=E n (r-r e , Q, 0-(H ra (r-r o , <) e N (r)> H (r) (8) 
is calculated (S 1 1 , S 1 2, S 1 3). hi the above 

e K(r) = e jll (r)xe„ 2 (r) (9) 

where die symbol "x" indicates a vector product. In this instance, an independent 
component in die projective component is one, diat is, oidy the magnitude of die 
vector. 

Now, let the coordinates of die magnetic field source be represented by a 
vector r c (x, y, z) and its angle of orientation by a vector 0^ 9 y , 9 2 ). In die 
following description, 9 X , 9 y and 9 Z will be referred to as an angle of rotation, a tilt 
angle and an azinurth, respectively. When die signal magnetic field is axiaUy 
symmetric, die axis of symmetry is regarded as die x- axis, and die angle of 
orientation is set to a vector 9 c (9 y , 9^). 

Further, according to die present invention, it is possible to calculate the 
magnitude H m (r-r c , 0 C ) of each component of the original signal magnetic field 
vector H m (r-r c , 0 C , t) by synchronous detection of the measured magnetic field 
vector H m (r-r c , 9 C , t) tiirough the use of a proper one component in Eq. (1) or (6) 
tiiat is a projective component. 

H.(r-r.. 0>(H m (r.r , t)Hj(r-r e , 0), (10) 

In die above, H m p (r-r c , 9 C , t) is a proper component of the projective 
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component given in Eq. (1) or (6). 

Fig. 4 shows the case where the number of noise magnetic field sources is 
essentially one. The processing of Fig. 4 is repeated at each measuring point. Li 
the case of essentially two noise magnetic fields, too, the processing shown in Fig. 

5 is repeated at each measuring point as is the case with Fig. 4. 

The magnitude H m (r-r c , 0 C ) of the original signal magnetic field H m (r-r c , 6 C , 
t) can be obtained by calculating Eq. (10) after calculating the projective 
component at each measuring point.. Concretely, die processing depicted in Fig. 

6 or 7 is carried out [(SI I, S12, S14, S15, S16) or (Sll, S12, S14a, SI 6)]. 

How to calcul ate the p ositi on and ang le of orientation: 
(1) Definitions of die coordina te system and t he angle of orientation- 
A definition will be given first of the coordinate system necessary for the 
description of the invention. 

A coordinate system is set which is fixed to the earth with die z- axis in a 
vertical direction (upward), which system will hereinafter be refeired to as a 
measuring coordinate system vector C M . The x- and y- axes are properly set so 
diat diey fomi a right-hand system. For example, they are set in parallel to die 
direction in which the side of a measuring frame is projected on a horizontal plane. 

This has for its object to calculate die coordinate vector r e (x, y, z) and 
angle-of orientation vector 9,(0* 9 y , 0 Z ) of the magnetic field source in die 
coordinate system. 

On the odier hand, as the coordinate system vector C c of the magnetic field 
source, a coordinate system widi die axial direction of the magnetic field source as 
die x axis is set so that the y- and z- axes are horizontal and vertical (upward), 
respectively, when the magnetic field source is placed horizontally. 

The angle-of-orientation vector 0 C of die magnetic field source is defined as 
the angle of rotation between die measuring coordinate system vector C M and die 
coordinate system vector C c as described below. In the first place, a coordinate 
system Co parallel to die measuring coordinate system vector C M is turned by an 
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azimuth angle 0 Z about toe z- axis (in either one of the coordinate system). This 
coordinate system will be referred to as a coordinate system vector C c , . Next, the 
coordinate system vector C cl is turned by a tilt angle 9 y about the y- axis of the 
coordinate system vector C c , itself. This coordinate system will be referred to as 
a coordinate system vector C c2 . Further, (lie coordinate system vector C a is 
turned by an angle of rotation 9z about the x- axis of the coordinate system vector 
C c2 itself. The angle-of-orientation vector 6 C is determined so that the resulting 
coordinate system becomes the coordinate system vector C c . 

(2) Explanation of independent measurands and unknowns 
When toe number of noise magnetic fields is essentially one, toe 
measurement of toe magnetic field at one place will provide two independent 
measurands. When the number of noise magnetic fields is essentially two, the 
measurement of toe magnetic field at one place will provide one independent 
measurand. Further, in case of obtaining components of toe original signal 
magnetic field by synchronous detection, the measurement of toe magnetic field at 
one place will provide three independent measurands. On the other hand, three 
coordinate components of the coordinate vector r c (x, y, z) of toe magnetic field 
source are unknowns. The azimuth angle 0 2 is not obtainable without a different 
criterion such as the direction of toe earth magnetism, and hence it is also an 
unknown unless toe digging head is provided with a direction sensor. When the 
digging head is made of a magnetic material, or when a buried steel pipe or similar 
magnetic body lies near the digging head, an accurate direction cannot be obtained 
even if toe direction sensor is provided; hence, toe azimuth angle Q z is an unknown 
in many cases. The tilt angle 9 y can easily be detected by a tilt angle sensor that 
detects toe vertical direction, and hence it is known in many cases. The same is 
true of the angle of rotation 0 X . hi particular, when the signal magnetic field is 
axially symmetric, if toe axis of symmetry is regarded as toe x- axis, the angle of 
rotation G x becomes meaningless and hence can be ignored. 

At any rate, measured magnetic fields need only to be obtained at different 
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positions so that independent measurands larger in number than unknowns. 
Arrangement of the measuring sy stem- 

For example, as shown in Fig. 1, magnetic fields are measured using a 
required number of three-axis magnetic sensors disposed on the ground so that 
5 dieir relative positions are known. Since the direction of the noise magnetic field 
changes for each position of measurement, it is necessary to perform for each 
position the step (S21) for detecting the position of the noise magnetic field and the 
step (S22) for calculating the projective component of the measured magnetic field 
as shown in Fig. 8 that is a flowchart of the procedure for position determination 
10 by die present invention, hi case of calculating die component of the original 
signal magnetic field by synchronous detection, die processing dierefor (in a 
specified step) needs to be earned out for each position. 
Flow of measurement processing: 

A detailed description will be given of the detection (S21) of the direction 
15 of die noise magnetic field, the calculation (S23) of the position of die signal 
magnetic field source dirough the use of the projective component and the 
calculation of die position f of die signal magnetic field through die use of the 
signal magnetic field component 
[EMBODIMENTS] 
20 Embodime nt in die case of virtually one noise magnetic field: 

hi case of using die projective component, when the number of unknowns 
is Nu(>l) and die number of noise magnetic fields is virtually one, die projective 
component vector H m p (r-r c , 6 C , t) expressed by Eq. (1) is calculated from magnetic 
vectors H m (r-r c , 9 C , t) measured at Nu/2 or more different positions, and die 
25 position vector r c and the angle-of-orientation vector 0 C are determined (S31) as 
depicted in Fig. 9 so tiiat a projective component vector at each magnetic field 
sensing point is essentially equal to die projective component vector B/(t-t c , G c , t) 
of a tiieoretically calculated magnetic field of die projective component at each 
position of measurement, by which it is possible to detect (S32) die position and 
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orientation of the magnetic field source. 

In case of using the magnitude of the signal magnetic field synchronously 
detected by Hie projective component, the magnetic field is measured at Ny/3 or 
more different positions. When the number Nu(>l) of unknowns is even, the 
5 number of unknowns and the number of independent measurands can be made 
equal to each other; if setting magnetic field vector H m (r-r c , 0 C , t) is measured at 
positions N m =Nu/2, a N L number of such equations as given below need only to be 
solved. 

K, P (r k -r c , 0 c ,O),-", < , P (« k -r c , <?)=0,* = 1 N m ,q = \, 2. (11) 

10 where q = 1, 2 and represents two directions parallel to the plane of projection but 
parallel to each other (S33). Accordingly, H m>q p (r-r C) © c , t) and H e /(r-r c , 0 C , t), 
where q = 1, 2, are the magnitudes of q-direction components of the measured 
magnetic field vector H m (r-r c , 8 C , t) and the theoretical calculated magnetic field 
vector H e (r-r c , 0 C , t), respectively. Let it be assumed diat the vector 0 C represents 
1 5 any one of 0 X , 0 y , 0 & (0 y , 0 2 ), (0 & 0 X ), (0^ 0 y ), (0^ 0y> md ^ md ^ me vector ^ 

represents any one of x, y, z, (y, z), (z, x), (x, y), (x, y, z) and <J>, where <|) represents 
an empty set. 

For example, when unknowns are the position vector r c (x, y, z) and 
azimuth angle 0 2 of the magnetic field source, the magnetic field is to be measured 
20 at two different positions and four equations such as given below are solved, by 
which it is possible to obtain (S34) the position vector r c (x, y, z) and azimuth angle 
0 Z of the magnetic field source. 

KVfo-";, 0 z ,O),-#, q P ('- k -r c , 0 Z ) = O, k = \, 2-q=\, 2. (12) 

hi the above, < . > t represents a time average. 
25 As depicted in Fig. 10, when the magnetic field is measured at N m (>Nu/2) 

different places more than Nu/2 in excess of die number Nu of unknowns, since a 
larger number of independent measurands than the number of unknowns can be 
obtained (S41, S42, S43), the position r c (e, y, z) and the angle of orientation 0 2 are 
calculated which provide 

15 



SS{||X fa /(* - '<> «• <D ( ~ /(ft " ft, * c | (13) 
where < > t represents the time average and min { .} means that vectors r c and 0 C 

are changed to obtain vectors r c and 9 C that minimize the contents in { . } . Further, 
a symbol w M indicates weighting. Hie Eq. (13) can be replaced with the 
following equations (S44). 




Assume diat vectors r c and 0 C in Eqs. (13), (14), (15) and (16) have the same 
meaning as in the case of Eq. (11). 

In the above description it does not matter whether the measured magnetic 
field vector H m (r-r c , G c , t) is a signal having passed through a band pass filter that 
permits the passage theredirough of only components close to die frequency of die 
signal magnetic field or a wide-band signal that is inhibited from die passage 
dirough die band pass filter, but die use of me signal having passed tiirougli die 
band pass filter increases die possibility of determining die position of die 
magnetic field with high reliability. 

(Embodiment in die case of essentially two noise magnetic fields) 
When die number of unknowns is Nu (>1) and the number of noise 
magnetic fields is virtually two, die projective component vector H m p (r-r c , 0 C , t) 
expressed by Eq. (4) is calculated from magnetic vectors H m (r-r c , G c , t) measured 
at Nu/2 or more different positions, and die position vector r c and the 
angle-of-orientation vector G c are determined (S51, S52) as depicted in Fig. 11 so 
that die above-mentioned projective component vector obtained at each position of 
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measurement essentially matches with the projective component vector H e p (r-r c , 0 C , 
t) of a theoretically calculated magnetic field at each position of measurement, by 
which it is possible to detect (S53) the position and orientation of the magnetic 
field source. 

5 hi diis case, the number Nu (>1) of unknowns and the number of 

independent measurands can be made to be equal to each other without fail; the 
magnetic field vector H m (r-r c , 9 C , t) is measured at different positions of die same 
number as diat N u number of unknowns, and Nu number of such equations given 
below need only to be solved. 

10 K P (r k -r CJ 0 C> /)), -///(r k -r CJ O c )=0, k = \ N u . (17) 

In this case-symbols H m p (r-r c , 9 C) t) and H e p (r-r c , 0 C , t) are the respective 
magnitudes of die projective component vector H m p (r-Tc, 9 C , t) and die vector 
H e p (r-r c , 9 C , t) of the measured magnetic field vector H m (r-r c , 9 C , t) and die 
dieoretical calculated magnetic field vector H e (r-r c , 9 C , t). Let it be assumed diat 
1 5 die vector 9 C represents any one of 9„ 0 y , 9 & (9 y , 9,), (9 a 9 X ), (9 X , 0 y ), (9,. 9 y , 9 Z ) 
and § and that the vector r c represents any one of x, y, z, (y, z) 5 (z, x), (x, y), (x, y, z) 
and <J), where <j> represents an empty set. 

For example, when unknowns are the position vector r c (x, y, z), azimudi 
angle 9 2 and tilt angle 9 y of die magnetic field source, the magnetic field is to be 
20 measured at five different positions and four such equations given below are 
solved, by which it is possible to obtain die position vector r c (x, y, z), azimutii 
angle 9 Z and tilt angle 9 y of the magnetic field source. 

(tf m P (r t - r c , O c , t)), - K P (r k - r C) ^ c )= 0, k = 1 5. (18) 

hi die above, < . > t represents a time average. And die vector 9 C = 9 c (9 y , 9J 
25 As depicted in Fig. 12, when die magnetic field is measured at N m (>Nu) 

different places more dian Nu in excess of the number Nu of unknowns, since a 
larger number of independent measurands dian the number of unknowns can be 
obtained (S61, S62), the position r c (e, y, z) and the angle of orientation 9 Z are 
calculated which provide 
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£^ P (r h -r..4.^ | -///(r k -r„^ 



(19) 



where < . > t represents die time average and min { .} means mat vectors r c and 9, 



r.,0, 



are changed to obtain r c and 0 C lliat minimize the contents in { . }. Further, a 
symbol w M indicates weighting. Hie Eq. (19) can be replaced with the following 
5 equations (S63). 

™j{|^- P (p k -r.,4, if) f - HJh - r „ <9 C |. (20) 

?^Jk P 0- k -r.. «. <J),-"/(r k -r c , 5>)f j (21) 
min[|w t |^/(r k -r c ^ c> ^ - ^ _ ^ ^ j (22) 

Assume tliatr c and 0 C in Eqs. (19), (20), (21) and (22) have the same meaning as in 
10 the case of Eq. (16). 

hi the above description, it does not matter whether the measured magnetic 
field vector H m (r-r c , 9 C , t) is a signal having passed through a band pass filter that 
permits die passage therethrough of only components close to the frequency of the 
signal magnetic field, or a wide-band signal that is inhibited from the passage 

15 through the band pass filter, but the use of die signal having passed through the 
band pass filter increases die possibility of detemiining die position of die 
magnetic field with high reliability. 

Further, in case of using the signal magnetic field component obtained by 
synchronous detection, die direction of die noise magnetic field is determined and 

20 is used to obtain die projective magnetic field; and the subsequent processing is 
common to die cases of one and two noise magnetic fields. The projective 
component vectors H m p (r-r c , 0 C , t) at respective places are calculated from 
magnetic field vectors H m (r-r c , 0 C , t) measured at Ny/3 or more different places, 
and a proper one of the calculated vectors is used as a reference signal to perform 
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synchronous detection of the measured magnetic field vectors H m (r-r c , 0 C , t), 
thereby obtaining the magnitude H s (r-r c , 0 C ) of die original signal magnetic field 
component. By determining die position vector r c and the angle-of-orientation 
vector G c of the magnetic field source so mat the magnitude of the original signal 
5 magnetic field component and the magnitude H c (r-r c , 0 C ) of Ihe theoretical signal 
component are equal to each other, it is possible to detect the position and 
orientation of the magnetic field source. 

When the number N 0 (>l) of unknowns is a multiple of 3, die number of 
unknowns and the number of independent measurands can be made to be equal to 
10 each other; if the magnetic field vector H m (r-r Cj 0 C , t) is measured at positions 
N m =N„/Nu/3, Nu number of such equations as given below need only to be solved. 

H s (r k -r c , <0-H 0 (r k -r e> 0 c )= 0 , k = 1, ... ,N n . (23) 
Fig. 13 shows the flow of processing (S61, S62, S65). Here, H s (r-r c , 0 C ) and 
Hc(r-r c , 0 C , t) are die signal magnetic field calculated from the measured magnetic 
15 field vectors H m (r-r c , 0 & t) and the magnitude of die theoretical calculated 
magnetic field vector, respectively The magnitude Hs(r-r c , 0 C ) of the signal 
magnetic field vector is an averaged quantity as already explained. Let it be 
assumed diat die vector 0 C represents any one of 0 to 0 y , 0 a (0 y , 0 2 ), (0,, 0^ (0^ e y ) s 
(Ox, 6 y , 0 2 ) and (|) and that the vector r c represents any one of x, y, z, (y, z), (z, x), (x, 
20 y), (x, y, z) and ty, where represents an empty set. 

For example, when an unknown is die position vector r c (x, y, z) of the 
magnetic field source, the magnetic field is to be measured at one place and three 
such equations as given below are solved, by which it is possible to obtain the 
position vector r c (x, y, z) and azimudi angle 0 Z of die magnetic field source. 

25 H,(r-,;)-H c (r-r c )=0. (24 ) 

The vector r is the vector of the position of measurement 

When die magnetic field is measured at N m (>Nu/2) different places more 
dian Nu in excess of the number Ny of unknowns, since a larger number of 
independent measurands than the number of unknowns can be obtained, die 
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position r c (e, y, z) and the angle of orientation 0 Z are calculated which provide 




where min { .} means that vectors r c and 0 C are changed to obtain vectors r c and 9 C 

that minimize die contents in { }. Further, a symbol w M indicates weighting. 
The Eq. (25) can be replaced with the following equations. 

™n{!"4 H »k - r c , ^l-|K(r k - r c , *J)j. (26) 
™*{%™*\H s (r k -r B , ^)-H c (r k -r c) *jf J. ( 27 ) 

Fig. 14 shows die flow of processing (S61, S64, S66). 

A description will be given below of how to determine the direction of die 
noise magnetic field in the two embodiments described above. 

(How to determine die direction of the noise magnetic field in the case of 
virtually one noise magnetic field) 

(First mejhod) 

The first metiiod for determining the direction of the noise magnetic field is 
a metiiod where, in case of no signal magnetic field, a noise magnetic field is 
measured tiirough me use of the same measuring system as that for measuring the 
signal magnetic field. This situation is such as shown in Fig. 15 in which the 
detection of die digging position is disturbed by a noise magnetic field source lying 
in die vicinity of die digging route. Anodier case is that where the signal 
magnetic field source is equipped widi a function of receiving a command sent, for 
example, from die ground by some means and responsive to die command to stop 
die generation of die signal magnetic field. 

In tiiis instance, as depicted in Fig. 16, letting the measured magnetic field 
be represented by a vector H m (r-r c , 0 C , t), since this is essentially a noise magnetic 
field vector H n (r, t) as shown at a step (S71), an average value of its absolute values 
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can be used (S72) to calculate a direction vector, e n (r)=e n , x (r), ev y (r), e^r)), of die 
noise magnetic field as follows: 

K«( r - r - *» 

Alternatively, a root-mean-square value of tlie above absolute values can be used to 
determine the direction of the noise magnetic field by 

V(fe-(r t -r.,g.,Q)') | 

In tiiis case, a symbol H m>a (r-r c , 0 C , t) is an a component of the measured magnetic 
field (the noise magnetic field), and a is any one of x, y and z. 

hi the above description it does not matter whether the measured magnetic 
field vector H ltl (r-r c , 0 C , t) is a signal having passed through a band pass filter that 
permits me passage therethrough of only components close to the frequency of the 
signal magnetic field or a wide-band signal that is inhibited from the passage 
through the band pass filter, but the use of the signal having passed through the 
band pass filter increases the possibility of detemiining the position of the 
magnetic field with high reliability. 

(Second mefliod^) 

Stepl 

As shown in Fig. 17, at first , the frequency spectrum H m (a>) of the 
measured magnetic field vector H m (r-r c , 9 C , t) is calculated (S81, S82, S83) by the 
following equation. 

H M) = F (H„,(r k - r 0 , /)) (30) 

In this case, a symbol F( . ) represents a Fourier transform; the three components x, 
y and z of the measured magnetic field vector HUr-rc, 8 C , t) are each Fourier 
transformed, hi practice, the above-mentioned frequency spectrum can be 
calculated by FFT (fast Fourier transform) or the like of sampled values of the 
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measured magnetic field vector H m (r-r c , 0 C) t). 
Step 2 

Next, an angular frequency co b where i = 1, 2, N s , of a large-amplitude 
component, such as a line spectrum, is selected (S84) from the absolute value 
5 |H m (G))| of the frequency spectrum. For the component of each angular frequency 
©i, where i = 1, 2, . . . , Nns , a candidate unit vector e„(r, a)j), where i = 1, 2, . , of 
the direction of die noise magnetic field, is calculated (S85) by the method (1) or 
(2) described below. 

(1) The absolute values of the Fourier-transfonned x, y, and z components 
1 0 of die angular frequency concerned are used to calculate the candidate unit vector 
e„(r, G>i), where i - 1» 2, . . . , N» of the direction of the noise magnetic field, by die 
following procedures: 

I a=x, y, z 

where H,^), where a = x, y, z and i = 1, 2, N., is a ©j component by the 
1 5 Fourier transform of the a component of die measured magnetic field. 

(2) A narrow-band filter is formed whose pass band uses, as die center 
frequency, the angular frequency where i = 1, 2, .... Nns , and the same method 
as by Eq. (21) or (22) is used to calculate the candidate unit vector e„(r, a>j), where i 
= 1, 2, .... N B . That is, the candidate unit vector e^r, ©;)= (e^r), e^r), e^r)), 
20 where i = 1, 2, . . . , H n ,is calculated by 



or by 



, a=x,y,z;i=l...,N Bs . (32) 



f *=x f y,x ' I 



a-jc.j,, z;/ = l jyr.. (33) 
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Step 3 

The candidate unit vector e n (r, o>i)= (e ivc (r), e^r), e^r)), where i = 1, 2, ... , 
, is considered as the direction vector e n (r) of tlie noise magnetic field, and for 
each angular frequency where i = 1, 2, ... , , the same method as by Eq. (1) 
is used to calculate a projective component vector H m p (r-r c , 0 C , ©j, 1). 
H m P (r-r c , 0 e , a l ,i)=Ujr-r.,O t , () 

~ ( H «,( r - , 0 C , /)-e.(r. fl»,)>.(r, <»,) i = I, ... , N m . (34) 
The reason for which the angular frequency G)j is contained as a variable of die 
projective component is to explicitly point out that the projective component is 
dependent on die angular frequency ©j, where i = 1, 2, N 11S . A roper time 
interval T tesU which consists of NU durations where k = 1, 2, N«, each 
having a short time length At, is chosen, and die variation of the projective 
component vector H m p (r-r c , 0 C> © b t) for each duration T^ where k = 1, 2, . . N lest , 
is evaluated (S86). Assume, here, that each duration where k = 1, 2, . . . , N tesl 
does not overlap otiier durations. Concretely, a variance of N test statistics v wal>k (fflj) 
of die NU, where k= I, .... N^, wliich are calculated by any one of die methods 
described below, is calculated. 

(1) One or both of die means of absolute values of two orthogonal components by 

v .™ ll c(® i )=K,q P ( r - r c» 3> 'D Tik .?=l 2;*=l,...,A^;i = l, ...,N m . (35) 

where Otu represents die mean value in die duration Ty, and is a statistic 
calculated for the duration T^. 

(2) Mean of absolute values 

k fe)=dH m P (r-r c> a> h ty^, k=], ... , N tcst ; i= 1, ... , N ns . (36) 

(3) One or bodi of die means of squares of two orthogonal components by 
v ^ l k)=((^ m ., P (r-r c , 0 e ,a.„t)j) ,q = 1, 2; k = 1, ... , N^; i= 1 

(37) 
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(4) One or both of square roots of the means of squares of two orthogonal 
components by 

v .«i kW«^(#* ,"(••-»■.. 3. *>i> Oj) • 9 = 1, 2; * = 1, ... , / = 1, ... , N as . 

(38) 

5 For the statistics v eval , k , where k = 1, .... 1SU, calculated by these equations, the 
following equation 



Jmean k (jv evill k (aj.)- mean k (v evDl , k (^))) 2 ) 
var(a):)= — 7 ^ — ^ /- i 7/ 



(39) 



is calculated to obtain (S86) a value of (o ltinin that is the angular frequency a)j which 
minimizes var(©j). In the above, mean k (.) indicates averaging for die suffix k, 
10 that is, 

£0 



mean 



*(.)-*tr-. (40) 



The magnetic field of tlie angular frequency derives from the noise 
magnetic field, and tire direction of tire noise magnetic field becomes a vector e„(r, 

Wi,mi„). 

15 Incidentally, die angular frequency ©j >min , which ininimizes var((Di), needs 

only to be measured at one place and need not be obtained at every place where to 
measure die magnetic field. 

With this method, it is also possible to calculate a fluctuation in the 
direction of a vector H m p (r-r c , 8 C , co ia t) as well as an amplitude fluctuation given by 
20 Eq. (39) and to select die angular frequency © iimin at which the direction fluctuation 
becomes minimum or smaller dian a predetermined value. 

Incidentally, in the above description, die measured magnetic field vector 
H m (r-r c , 9 C , t) in Step 1 is a wide-band signal, and in Step 3 it does not matter 
whether the measured magnetic field vector H m (r-r c , 9 C , t) is a signal having passed 
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through a band pass filter that permits the passage therethrough of only 
components close to the frequency of the signal magnetic field or a wide-band 
signal diat is inhibited from the passage through die band pass filter, but the use of 
the signal having passed through the band pass filter increases the possibility of 
5 detennining the position of die magnetic field with high reliability. 

Fig. 18 shows how to select frequencies fi(=a>,/27i), f 2 (=a> 2 /27i), 
fn(=(<V27i, where n=N ns ), at which the frequency spectrum becomes maximum; 
Fig. 19 shows die flow of processing including step S91 for obtaining the 
candidate vector; and Figs. 20 and 21 show die flow of processing including steps 
10 S101 and S102 or steps SIM and SI 12 for evaluating the candidate vector and for 
detecting die direction of die noise magnetic field. 

The frequency spectrum H m (o)) need not always be used. That is, die 
signal magnetic field is periodically turned OFF/ON following a predetemiined 
procedure; die period T^od is divided into equally-spaced durations; die candidate 
1 5 unit vector e n (r, tj), where i = 1, . . . , N re , is used in place of die candidate unit vector 
e n (r, ©j), where i = 1, N„s, which is calculated by Eqs. (25), (26) and (27); and 
diereafter, the duration diat minimized the variance by Eq. (33) is calculated by die 
processing described above. By tiiis, die vector e,,(r, tj) in that duration can be 
adopted as die direction of he noise magnetic field. 
20 (Third mediodl 

hi the second mediod, the large-amplitude angular frequency ©;, where i = 
1, 2, N s , is selected from die absolute value |H m (co)| of me frequency spectrum 
H m (a)), but the vector e n (r, Q) i>min ) can be obtained as the direction of die noise 
magnetic field in exactly (lie same manner as in die case of die second mediod, by 
25 selecting a proper frequency band neighboring die frequency of die signal 
magnetic field, setting properly-spaced test frequencies free from the frequency of 
the signal magnetic field in the frequency band and regarding die test frequencies 
as die angular frequency ©j in die second mediod. 

As is die case widi die second mediod, die angular frequency © ijmin at 
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which var(a)j) becomes minimum needs only to be obtained. 

Fig. 22 shows how to select the frequencies f,(=®,/2re), t 2 (=(D 2 /2n), .... 
ttKNA, where n=N ns ) at which tlie frequency spectrum becomes maximum! 
The flow of the subsequent processing is the same as depicted in Figs 17 19 and 
20. 

(Fourth method) 

Tlie signal magnetic field is periodically stopped under control of a 
predetermined procedure. For example, tlie signal magnetic field is periodically 
stopped by a predetermined time interval. Since the intensity of tlie magnetic 
field being measured decreases while tlie signal magnetic field is stopped, the OFF 
period of tlie signal magnetic field is identified by regarding the 
intensity-decreasing period essentially as (lie predetermined OFF period, and the 
direction of tlie magnetic field measured during (he OFF period is used as tlie 
direction of the noise magnetic field. The direction of the noise magnetic field 
can be obtained using tlie same method as tlie first one. Fig. 23 shows how tlie 
amplitude of die measured magnetic field in this method varies with time. 
(Fifth method) 

Tlie signal magnetic field is periodically stopped following a 
predetermined procedure. This is carried out as described below under (1) and 
20 (2). 

(1) To stop the signal magnetic field on a rectangular-wave-wise: 

Tlie signal magnetic field is repeatedly turned ON and OFF witii the period 

Tperiocb 

for instance. 
s(t) = l, i k zt<t k +t aap . 

= - 1 > ( 41 ) 
and when a sequence s(t) is 1, die signal magnetic field is turned OFF, but when 
the sequence is -1 , the signal magnetic field is turned ON. In the above, 
'k+i -'k = period, k= 1, 2, 3, ... ^42) 
(2) To stop die signal magnetic field on a pseudo-random-signal- 
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-wise 
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basis: 

For example, when the value is "-1" in a random sequence like an 
M-sequence consisting of unit periods T mlit of the same length Nm, die signal 
magnetic field is turned ON, but when the value is "1," the signal magnetic field is 
5 turned OFF; and this sequence is repeated. In this case, the time average of the 
sequence is set to 0. 

Fig. 24 shows temporal variations in the amplitude of the measured 
magnetic field in the case of the signal magnetic field being stopped by the method 
(1). Fig. 25 shows temporal variations in the amplitude of the measured magnetic 
10 field in die case of the signal magnetic field being stopped by the method (2). 

Next, as depicted in Fig. 26, die correlation function between the sequence 
s(t) and die norm of the measured magnetic field or die absolute value of its 
particular component is calculated (S121, S122, S123). As die correlation 
function, it is possible to use any one of diose calculated by 

15 /?(r)=f;; +//T ^|H m (r-r c> 0 c , $Q-r)dl. (43) 

^W=f^ TrpOW K«(r-r C5 0 C , /^-r)/U= W . (45) 

K(r)= ft*** 9 ** 3> '))> ' *>*• « = *.>v- (46) 

hi diis case, die period for which to detect die correlation is set to an integral 

20 multiple N T T pcr jod of die period T period . 

Hie OFF state of die signal magnetic field can be detected (S125) from the 
time T^ync at which any one of die above correlation functions become maximum 
(S124). That is, a sequence s^-t^), which starts at the time x=t sync , is used, and 
when die sequence s(t) is "1," die signal magnetic field is regarded as being OFF; 

25 in diis way, the ON/OFF operation of die signal magnetic field is determined. 

27 



In the thus determined signal magnetic field OFF period the direction of 
the measured magnetic field is detected, and die direction is regarded as the 
direction of the noise magnetic field (S126). When the sequence s(t) is "-1" (a 
second numerical value), the signal magnetic field is turned ON, and when the 
5 sequence is "1" (a first numerical value), the signal magnetic field is turned OFF; 
this sequence is repeated in this way. 

(Sixth method) 

As is the case with the fifth method, the signal magnetic field is 
periodically stopped under control of such a predetermined procedure as described 
10 below. 

(1) To stop the signal magnetic field on a rectangular-wave-wise: 

The signal magnetic field is repeatedly turned ON and OFF with the 
period for instance. 

*(/) = l. i k <t<l k +t sinp . 

(41) 

= -1, ',+/ sl0 p^'<W 
15 and when a sequence s(t) is 1 5 the signal magnetic field is turned OFF, but when 
die sequence is - 1 , the signal magnetic field is turned ON. In the above, 

"Ac = ^period* £ = 1> 2, 3, ... (42) 

(2) To stop the signal magnetic field on a pseudo-random-signal-wise 

basis: 

20 For example, when the value is in a random sequence like an 

M-sequence consisting of unit periods of the same length N M , the signal 
magnetic field is turned ON, but when the value is "1," the signal magnetic field is 
turned OFF; and this sequence is repeated. 

hi this case, the sequence s(t) is chosen so that it changes for each 
25 predetennined time unit At m n. And, the time average of the sequence is "0." 

Next, as depicted in Fig. 27, the correlation function between the sequence 
s(t) and the norm of the measured magnetid field or the absolute value of its 
particular component is calculated (S131, S132, S133) as in case of the fifth 
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method. As the correlation function, it is possible to use any one of those 
calculated by 



^)=C +WT7, > m (r-.-^ c ,4^-^- (43) 



R ^ = J, k V ( H ™ ( r - r - ')) - *■>*. 

5 «aW=f + " TU K <I (.-r tl ft.'J<'-r>*,a = w . 



(44) 
(45) 
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'4(0- £ ^^^^(r- Z O, t)J<(-r)d(, cc=x,y,z. (46) 

Iii diis case, the period for which to detect the correlation is set to an integral 
multiple NxTpeHod of die period T^. 

Li this Distance, diere are present, in general, plural times x=t wic j c (k = 1, 2, 
10 . . . , N sync ) in which the correlation function becomes maximum and the maximum 
value exceeds a predetermined value (S 134). Assume, for example, that (k 
= 1, 2, N^) is an arrangement of such times in order of time. When the 
correlation value between the sequence s(t) and the signal magnetic field is 
appropriate, 

A'sync k = 'sync k _ 'sync. 1 . * = 2 N^. (47) 

is virtually an integral multiple of die time unit AW Then, the average of the 
value resulting from die subtraction of an integral multiple M^AW of the time 
unit AL^, where k = 2, N^, from AW, where k = 2, .... N^, is 
calculated by 

S( A W-My„ckAC,) 

20 />/ - *° 2 

# ri • (48) 

lii this case, 

*sync = ''sync, 1 + ^syuc (49) 

provides die beginning of die sequence signal corresponding to the ON/OFF 
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operation of the signal magnetic field. 

Accordingly, the period during which the signal magnetic field is OFF can 
easily be set based on the sequence sO-t^). 

By applying the same method as the first method to die magnetic field 
vector H m (r-r c , 0 C , t) measured in this period, the direction e n (r) of the noise 
magnetic field vector H n (r, t) can be calculated (S138). 
(Seventh method) 

This mediod will be described below with reference to Figs. 28, 29, 30 and 

31. 

As is die case with die fifth mediod, die signal magnetic field is 
periodically stopped, for example, by such a procedure as described below 
(1) To stop die signal magnetic field on a rectangular-wave-wise: 
The signal magnetic field is repeatedly turned ON and OFF with die 
period T^od, for instance. 
15 <0 = l t*zi<t k +/, lop . 

and when a sequence s(t) is 1, the signal magnetic field is turned OFF, but when 
die sequence is -1, die signal magnetic field is turned ON. In the above, 
'k+i -'k - ZJeriod, k=\,2, 3, ... ^42) 
(2) To stop die signal magnetic field on a pseudo-random-signal-wise 

20 basis: 

For example, when die value is "-1" in a random sequence like an 
M-sequence consisting of unit periods T^ of die same length N M , the signal 
magnetic field is turned ON, but when the value is "1," die signal magnetic field is 
turned OFF; and diis sequence is repeated accordingly. In diis case, the time 
25 average of die sequence is "0." 

Next, die correlation function between die sequence s(t) and me measured 
magnetic field H m (r-r c , G c , t) is calculated (S141, S142, S143). The period T^ 
is divided into equally spaced N*, sections of a length T^, and either one ofthe 
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following calculations is conducted. 

; ^('*)= I'^K «('• " r c , 4(/ - t k )dl, k = 1, ... , W div ; a = x,y,z. (50) 

^('ic) = f ^VC^ «( r - r «- ^ 0) 2 4 - f k >*. *- 1, ... , ^ div ; a = (51) 

In tliis case, a symbol R^), where a = x, y, z and k = 1, N^, is the time 
5 correlation between an a component H nvx (r-r c , 9 C , t) of the measured magnetic 
field H m (r -r Cj 9 C , t) and Hie sequence s(t). And 

tt-to+k-T*. *-] N^. (52) 

The measured magnetic field may also be correlated with the time that is an m- 
multiple of the period T period . That is, 

10 /4(^)={ k+m7poind K a (r-i;, 0 C , t]s mp (t-l k )dl, k=l ...,N div ; a = x ,y,z. (53) 
R A)=\l +mT ^-i( H ^ r -^ 0)X('-'tM k = l, ...,N, iv ; a=x,y,z. (54) 

hi this instance, the sequence s(t) is replaced widi S^t), where 
5^(0 = 4+^) 

S mp (i) = .<<), 0<t<7^, rf . ( 55 ) 

This is followed by calculating the component H ra p (r-r c , 0 C , ^ t), where k = 
15 1, . . . , Ndj v , of the measured magnetic field H m (r-r c , 0 C , t) projected on the vector 
e„(tk)=(Rx(tk)> R y (t k ), R 2 (t k )), where k = 1, . . . , 
formed by correlation functions RJ&), where a = x, y, z , corresponding to the 
respective components x, y and z of die measured magnetic field H m (r-r c , 0 C , t). 
Here, a symbol t k contained as a variable of die projective component indicates that 
20 the projective component depends on the variable k 

hi tiiis mediod, die vector e^t*) can be used as the direction of the noise 
magnetic field through utilization of die tune t k in which a fluctuation in the 
absolute value of the projective component vector H m (r-r c , 9 C , ^ t), where k = 1, 
• • • 9 Ndjvj 
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. ) J((«- , fr-1-«.'..0-<|K.'fr-r..«./ t . / ^y 

(|H."(r-^, «,„,.<$, — = ' (56) 

becomes minimum or smaller than a predetemnned value (S145a). In the above 
(.), means die calculation of the time average 

5 Fu, * er - vm ^°f^x-c < , ra p„„e„ t H ra ,( Mc ,9 c ,t l ,0,y-co m po n e„t 
IMr-re, e„ (. t) and Z -compo„ent H,„,,(r-r 0 , 8, K t) of the projective component 
vector H„,(r-r c , 6 C , li, t), where k = 1 N fc , 

<l".'(--r,. (57) 
a- J :, X2 ;i.i, 

are calculated, and me veclor ^ can he used as the direction of the noise 
10 magnetic field through utilization of the time k in which the sun, of the 
above-mentioned variances 

.-?,.. V "" ( '' ) (58) 

or 



Jew 



(59) 



15 becomes minimum or smaller than a predetennined value 05145b) In this case 
.here is fce possibility that the con-elation functions RA ), ^ at ^ 
certatn tune k have lost their original signs. Hence, it is necessary to evaluate me 
fluctuation of fce projective component at each time for four combinations 

RAX W *A» [RA), R,( tt - Rl(tk)J , _ w> RA)] _ 
20 -R^lkX -Rrfyj. 

Further, the period during which the signal magnetic field is OFF can 
easdy be set (SI45c, S145d) baaen on the sequence s(t- tl ). By app.ying same 
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method as the first one to the magnetic field vector H n (r-r c , 0 C , t) measured in this 
period, the direction e„(r ) of the noise magnetic field vector H n (r, t) can be 
obtained (S 146). In the above, use is made of the time t* of the equally spaced 
Ndiv sections of a die time lengdi I* divided from the period but it is also 
5 possible to use the time when the correlation function given by Eq. (45) or (46) 
becomes maximum or when the correlation function becomes maximum and 
exceeds a predetermined value. 

(Mefliod fo r determining die direction of die noise magnetic field when the 
number of nois e magnetic fields is virtually two) 
10 Even in a case of two noise magnetic fields, when a first one of the two 

noise magnetic fields has far higher an intensity than the second noise magnetic 
field at a first frequency and the second noise magnetic field lias far higher an 
intensity than the first noise magnetic field at the second frequency, the directions 
of die fust and second noise magnetic field can easily be calculated dirough 
15 utilization of diese frequency components in die measured magnetic field. 

If the first or second frequency is close to the frequency of the signal 
magnetic field, die directions of the noise magnetic fields can be determined from 
the measured magnetic field having passed tiirough a band pass filter that permits 
the passage therethrough of only frequencies near tiiose of the signal magnetic field, 
20 by the same mediod as die fourth or fifth mediod for use in the case of virtually one 
noise magnetic field 

When eitiier of die first and second frequencies does not equal to the 
frequency of die signal magnetic field, die directions of the respective noise 
magnetic fields need only to be calculated by the same method as the first one for 
25 use in the case of virtually one noise magnetic field. 
(Odier embodiments) 

The present invention is also effective when the signal magnetic field 
generated by the magnetic field source is virtually axially symmetric, and the 
invention permits position determination with a smaller number of magnetic 
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sensors or by magnetic field sensing at a smaller number of positions than in the 
case of a magnetic field of low symmetry. 

Further, when only one noise magnetic field affects the measurement of die 
digging position and the tilt angle of die magnetic field source is known which is 
5 an inclination of die axis of symmetry corresponding to the axial direction of die 
signal magnetic field set in die magnetic field source widi respect to die vertical 
direction, die projective component of die magnetic field, measured at each of two 
or more different positions, on a plane perpendicular to the direction of die noise 
magnetic field sensed at each of the magnetic field sensing positions is calculated; 
1 0 the position of die magnetic field source and its azimuth angle that is die direction 
of the axis of symmetry in a horizontal plane can be calculated from die above-said 
projective component. 

Further, according to this invention method, when virtually one noise 
magnetic field affects the measurement of the digging position, the projective 
15 component of die magnetic field, measured at each of diree or more different 
positions, on a plane perpendicular to die direction of the noise magnetic field 
sensed at each of die magnetic field sensing positions is calculated; die position of 
the magnetic field source, its tilt angle tiiat is an inclination of die axis of symmetry 
corresponding to die axial direction of die signal magnetic field set in die magnetic 
20 field source with respect to the vertical direction, and die azimuth angle of die 
magnetic field source that is he direction of die axis of symmetry in the horizontal 
plane can be calculated from the above-said projective component. 

Further, when virtually two noise magnetic fields alone affect die 
measurement of die digging position and die tilt angle of the magnetic field source 
25 is known which is an inclination of the axis of symmetry corresponding to me axial 
direction of die signal magnetic field set in the magnetic field source witii respect 
to die vertical direction, die projective component of die magnetic field, measured 
at each of four or more different positions, on a straight line perpendicular to both 
of the direction of a first one of the two noise magnetic fields sensed at each 
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magnetic field sensing position and the direction of the remaining second noise 
magnetic field sensed at the same position is calculated; the position of the 
magnetic field source and its azimuth angle Iliat is the direction of the axis of 
symmetry in the horizontal plane can be calculated from die above-said projective 
5 component 

Further, when virtually two noise magnetic fields alone affect the 
measurement of die digging position, the projective component of the magnetic 
field, measured at each of five or more different positions, on a straight line 
perpendicular to both of the direction of a first one of die two noise magnetic fields 
10 sensed at each magnetic field sensing position and the direction of the remaining 
second noise magnetic field sensed at the same position is calculated; die position 
of the magnetic field source, its tilt angle fliat is an inclination of the axis of 
symmetry corresponding to the axial direction of the signal magnetic field set in 
die magnetic field source with respect to the vertical direction, and the azimuth 
15 angle that is die direction of die axis of symmetry in die horizontal plane can be 
calculated from the above-said projective component. 

Moreover, when virtually two noise magnetic fields alone affect the 
measurement of die digging position, die frequency component of a first one of die 
two noise magnetic fields, in the vicinity of which die remaining second noise 
20 magnetic field aid die signal magnetic field have substantially no frequency 
components, is measured to diereby permit detection of the direction of die first 
noise magnetic field in tenns of vector; and die frequency component of the 
second noise magnetic fields, in die vicinity of which the first noise magnetic field 
and die signal magnetic field have substantially no frequency components, is 
25 measured to thereby permit detection of die direction of die second noise magnetic 
field in tenns of vector 

In die present invention, it is effective to use, as a magnetic sensor, a 
three-axis magnetic sensor that senses three magnetic fields orthogonal to one 
anodier at substantially the same position. 
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The magnetic sensor for use in die present invention may be any lands of 
sensors as long as they are capable to three magnetic fields orthogonal to one 
another at substantially the same position, but the three-axis magnetic sensor is 
suitable which senses three magnetic fields orthogonal to one another at 
5 substantially the same position. Alternatively, it is possible that one magnetic 
sensor capable of sensing a magnetic field in only one direction is turned at the 
same position toward three orthogonal directions one after another to sense the 
three orthogonal magnetic fields. 

hi carrying out the present invention, it is possible to employ such a frame 

10 12 as depicted in Fig. 32 which has magnetic sensor fixing means 11 mounted 
thereon to fix three-axis magnetic sensors and a tilt angle gauge 13 for detecting 
the inclination of the frame with respect to a vertical direction. The position of 
each magnetic sensor fixing means on the frame is known, and the magnetic sensor 
fixing means possesses a function of fixing the magnetic sensor in a predetermined 

1 5 orientation to the frame. The magnetic sensor fixing means 11 is provided with, 
for example, three faces orthogonal to one another, and has a mechanism that fixes 
the magnetic sensor at a predetemiined angle when a predetemiined face of the 
sensor case is pressed against any one of the three faces of the sensor fixing means. 
One magnetic sensor or one three-axis magnetic sensor is fixed to these magnetic 

20 sensor fixing means one after another to sense Hie magnetic fields. 

hi another alternative, a plurality of magnetic sensors may be fixed in 
predetemiined orientations to the frame 12 at a plurality of positions to 
simultaneously sense magnetic fields at the plurality of positions. 

As described above, the present invention employs a frame provided with a 

25 plurality of magnetic sensor fixing means each capable of removably or fixedly 
mounting a three-axis magnetic sensor and a tilt angle sensor capable of detecting 
the tilt angle of an orthogonal coordinate system to Hie vertical direction, the 
magnetic fixing means being mounted on Hie frame so that their positions and 
orientations are known; (lie magnetic sensor is removably or fixedly mounted on a 
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required number of magnetic sensor fixing means to sense magnetic fields, and the 
tilt angle of the frame during magnetic field sensing and the orientation of the 
magnetic sensor at each magnetic sensor mounting position with respect to the 
frame are used to calculate from the magnetic field sensed at each magnetic sensor 
mounting position the sensed magnetic field, a noise magnetic field and a signal 
magnetic field as vectors in a coordinate system fixed to the ground. 

The magnetic field generating means for the signal magnetic field in the 
present invention may be a coil. The magnetic field generating means may be 
one electric wire as well, or may also be one electric wire that is straight only in the 
vicinity of the place of position determination. 

As described above, even if a buried power line, railroad tracks, or similar 
noise magnetic sources are present near a construction site, the present invention 
permits highly reliable position determination without being affected by noise 
magnetic fields generated by such noise magnetic field sources. 

The present invention is intended for measuring die digging position in the 
non-open-cut method of excavation, but is applicable as well to many technical 
fields that involve position determination by sensing magnetic fields. 
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